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Abstract 

We study two pion correlations from Skyrmion and antiSkyrmion collision, us¬ 
ing the product ansatz and an approximate random grooming method for nu¬ 
cleon projection. The spatial-isospin coupling inherent in the Skyrrne model, 
along with empirical averages, leads to correlations not only among pions of 
like charges but also among unlike charge types. 
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I. INTRODUCTION 


We have studied nucleon-antinucleon annihilation at rest in the large N c limit of QCD 
in a series of papers OKI- The gross features of annihilation emerge correctly from this 
picture. These features include the pion momentum distribution, number and charge spec¬ 
tra and branching ratios into meson channels. In these studies the initial baryon density 
distribution was taken as spherically symmetric and only the subsequent evolution of the 
mesons was treated dynamically. In the large N c limit this dynamics is the classical dynam¬ 
ics of Skyrrne like models. To discuss annihilation in flight, pion correlations, the effects 
of the spin and isospin of the individual baryons and the initial state dynamics involves 
considerable improvement and sophistication of our previous starting point. We also need 
to address the question of quantum, or finite Nq, corrections. These are usually introduced 
by rotating the Hedgehog (and antiHedgehog) slowly and quantizing the rotational degrees 
of freedom. In the Skyrrne picture of NN collision, we need to collide a rotating Skyrmion 
with a rotating antiSkyrmion, project out the nucleon and antinucleon states, and follow the 
dynamical evolution of the classical system to pion radiation. This is, of course, a much more 
demanding task than the study of classical dynamics of annihilation without grooming [f||, 
which is already an elaborate and numerically challenging calculation. These difficult issues 
of classical Skyrmion dynamics and quantum correction in annihilation bear strongly on 
pion Bose-Einstein correlations ||, charge asymmetry ||.|7[] in pp —> 7r + 7r“ and annihilation 
in flight in general. 

This paper is a first modest step in the study of annihilation in flight and of pion intensity 
correlations. We study annihilation in flight with the product ansatz and with no subsequent 
pion dynamics. Pion intensity correlations emerge as a consequence of random grooming of 
approaching Skyrmion and antiSkyrmion. The result are similar to the commonly considered 
Bose-Einstein effects resulting from thermal decoherence. However, the correlation in the 
present study results from the spatial-isospin coupling inherent in Skyrrne soliton. Pion 
correlations from this spatial-isospin coupling were studied by Blaizot and Diakonov Q in 
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heavy ion collision. 

In the next section, we obtain slowly moving Skyrmion configurations from Galilean 
transformation. We approximate the full dynamics by the configuration of a Skyrmion and 
antiSkyrmion nearing each other. Such a configuration at the point of closest approach is 
used in Section 3 to obtain a coherent state amplitude describing the pion radiation from 
annihilation. Interaction between the Skyrmion and antiSkyrmion or among the pions is not 
considered. In section 4, we discuss the random grooming of the Skyrmion and antiSkyrmion 
sources and resulting pion correlations. Section 5 outlines the case of Lorentz transformed 
Skyrmions and compares the results with those from the Galilean case. Section 6 gives our 
results and prospect for future development. The various approximations we make, some of 
them quite drastic, are discussed in the sections where they are employed. 

II. CLASSICAL MOVING HEDGEHOGS: GALILEAN CASE 

The large N c dynamics we consider is the Skyrrne model |9j which is a classical nonlin¬ 
ear field theory of pions only in which the baryons emerge as nonperturbative topological 
solitons. These baryons are not nucleons but rather hedgehogs (so called because the pion 
isospin and spatial directions are correlated). The hedgehogs are a superposition of all possi¬ 
ble baryonic states with isospin equal to spin. The nucleon, or antinucleon, is then projected 
out from the hedgehog configuration by grooming or rotating the hedgehog. When a hedge¬ 
hog and antihedgehog approach each other, there is a complicated interaction coming from 
the nonlinear nature of the Skyrrne dynamics, leading finally to pion radiation. This at least 
is what happens for close collisions, and has been seen in the one full calculation that has 
been done. [jj The nonlinear nature of the Skyrrne lagrangian and the presence of fourth 
order derivatives makes the calculation numerically unstable. Thus a full investigation of the 
nature of Skyrmion antiSkyrmion collisions as a function of energy and impact parameter 
and of the subsequent pion radiation has not yet been carried out even for the hedgehog 
anti-hedgehog orientation, let alone for groomed hedgehogs. We do not even try. Rather 
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we take as the initial state a hedgehog and antihedgehog moving in opposite directions and 
at their point of closest approach. This is the product ansatz for moving hedgehogs. The 
hedgehogs themselves we take as solutions of the Skyrrne equations. We then use their 
classical pion configuration as the source of pions for our quantum coherent state with no 
further dynamics among those pions. The calculation of the Cal Tech group Pf and the 


one model calculation we carried out in simple geometry |T0| suggest that once annihilation 
begins, it proceeds very rapidly, at the causal limit, and that the subsequent pion wave in 
the radiation zone looks very much like the wave one would obtain from a free wave fourier 
transform of the source. We use these facts and the fact that anything else is too difficult 
for now to justify our simple approximations. 

In the center of mass, we represent the colliding nucleon and antinucleon by a classical 
Skyrmion and antiSkyrmion approaching each other with velocity v and impact parameter 

We make the approximation that annihilation and the subsequent pion radiation comes 
at the instance of closest approach, as shown in Fig. 1. Final state interactions among the 
resulting pions are neglected as discussed above. 

The pion held configuration of a moving hedgehog is the Lorentz boosted static hedgehog. 
For small velocities, we can approximate it by a Galilei boosted hedgehog, which we study 
first. 

The classical held of pions in Fig. 1 can be written as 
<3>(r,t) = $ a (r, t) + 

<F a (r, t) — (r - vt + b/2)g(\r - vt + b/ 2|) 

<f> fe (r, t) = —(r + vt - b/2)g(\r + vt-b/2\) (1) 

where g(r) = f(r)/r (/(r) is the usual Skyrrne prohle function, modulo /„■), a is the Skyrmion 
and b is the antiSkyrmion. The fact that the iso-vector holds on the left in (!]) are written 
in terms of radial spatial vectors on the right is the hedgehog feature. 


Bold face letters denote 3-vectors. Overhead arrows denote iso-vectors. 
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III. REQUANTIZING 


To form a coherent state, we need to decompose the classical pion field ([I]) into a plane 
wave and pick out the positive frequency part. We neglect the pion interactions and thus 
the pion wave satisfies the free wave equation (the isospin index is suppressed for now), 

VV(r, t) = (+ m 2 )^(r, t) (2) 

where m is the pion mass. 

We introduce the momentum space amplitude h(k, u) by 

ip(r,t) = J d 3 kdue lk ' r e~ lujt h(k,uj) (3) 

and we find 


h(k, u) = 5(k~ + m — u )G(k, u) 


(4) 


and therefore 


i/>(r,t) = f d 3 kdu(e ik ■ r e- i “ kt G{k,u Jk )/(2uj k ) + e~ tk r e iuJkt G(k, -uj k )/(2uj k )). (5) 


The positive frequency part of the wave comes from the first term, the one with e lUkt and 


therefore we need G(k, u k ), [oj k = Vk 2 + m 2 ). From the field and its rate of change at t — 0, 
the instant of nearest approach, we obtain 


^(r*, 0) = / d 3 ke 


3i „ik-r G(k,u> k ) G(k, uj k ) 


2(jO k 


= J d 3 ke‘ kr 't I (k) 


( 6 ) 


and 


dtp(r, 0) = r d 3 hc ik-r ^ „■) ^k) - G(k, -u k ) 


dt 


= J d 3 ke lk ' r ^ 2 (k) 


(7) 


We therefore see that 
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G(k, uj k ) = Uk'yi(k) + * 72 (fe) 


(8) 


and the fourier amplitude of the held that we need is given by 

G{k,u k ) 

h{k,u k ) = —-- (9) 

Zujk 

To obtain 71 , 72 and thus G(k), we fourier transform the held conhguration (jl]) and its 
time derivative 


7i(fc) = 


(2vr) 


<Pre- ik ' r $(r, Q) 


( 10 ) 


72 (k) = 


d 3 re -ikrdi { rQ) 


( 11 ) 


(27t) 3 J dt 

We split the amplitude into its contributions from the hedgehog and from the antihedgehog. 


7 —> la + lb (12) 

*“=( 2 ^/*e-‘ fe > + |)s(|r + ||) (13) 

= “(2 \fj d3re ~‘ kT ^ - T 2 M\r - j|) (14) 

Introducing 

G(k 2 ) = j± j - 3 fd 3 pe-‘ k Pg(\p\) (15) 

it is easy to show that 

7 ia = 2ike lk ^G'(k 2 ) (16) 

7 i b = -2ike~ ik ^G\k 2 ) (17) 

with G'(k 2 ) = dG/d(k 2 ). 


$ a (r, 0 ) = —vg(\r + ^ | 2 ) - (r + |) « • 2 (r + s'(|r + ^ j 2 ) 
$&(**, 0 ) = -u^(|r - ^| 2 ) - (r - ^) v ■ 2{r - |) g\\r - ^| 2 ) 
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(18) 

(19) 








where 


ftp 2 ) = 


dg_ 

dp 2 


( 20 ) 


It is straightforward to show 


72 a 


1 

(2vr) 3 


d 3 re~ ik - r $ a (r, 0) 



G(fc 2 ) + ft 


( 21 ) 


where 

ft = d 3 re~' k r (r + \)[-2v ■ (r + \))g'{\r + *|) 

= i J d 3 pe~ ik 'Pp(-2v ■ p)g\p 2 ) 

= / d 3 pe- ik f) v ■ pS7 P g{p 2 ) 

= J d 3 pg(p 2 )V P [e- ik P v • p] 

= e ik ivG(k 2 ) + 2e ik ^k v ■ kG\k 2 ) ( 22 ) 


Note we use Vpp(p 2 ) = 2 pdj^ and then integrate by parts. Thus the fourier amplitude (^) 
of the held of the Skyrmion we need is 

ha(k, u k ) = ikft k G\k 2 )(l + ^-) (23) 

Similarly we have for the held from the antiSkyrmion 

h b (k, uj k ) = —ike~ lk ^G\k 2 )( 1 — — — ) (24) 

A quantum coherent state represent the configuration in Fig. 1 can then be written as 

\h) — exp[ J d 3 k(h a (k,u k )+ h b (k,u k )) ■ ad~(k)]\0) (25) 
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IV. RANDOM GROOMING AND CORRELATIONS 


We next consider grooming of the hedgehog and anti-hedgehog to introduce spin and 
isospin. In fact we will argue that we do not have to quantize the grooming but that we can 
consider random grooming. This is justified by consideration of time scales. Spin and isospin 
are introduced by spinning or rotating the hedgehog at a fixed rate. That rate of rotation 
is slow compared to the annihilation time. Thus during any given annihilation the groomed 
hedgehogs do not rotate much. But if there is no polarization, their orientation is random 
at the moment of annihilation and it is that random orientation or random grooming that 
we average over when we consider many events. 

As shown in the calculation of Skyrmion and antiSkyrmion annihilation PlJTOl, the pions 
emerge promptly from the interaction region. This happens in about the time for light 
signal to travel across the Skyrmion, which is T c — 1/m. Here m is the pion mass. The 
rotation period can be estimated from the N — A splitting which is of the order of m. The 
rotation energy is E r « Ir 2 oo 2 se M/(m 2 T 2 ) ~ m, thus T r = 1 /m^M/m. M is the nucleon 
mass. The ratio of rotation to annihilation time is then T r /T c ~ yjM/m ~ 3. If we were to 
estimate the rotation time by the angular momentum (of order 1, noting h — 1), we would 
have J ~ Iu> ~ M/m 2 /T r ~ 1 and T r = M/m 2 . The ratio would be T r /T c ~ M/m ~ 7. 
Both these estimates suggest that the approximation of taking a fixed orientation during 
the fast annihilation is a reasonable one. This is something that should be improved in more 
detailed calculations. In the limit of slow rotation, the isospin orientation of each Skyrmion 
can be considered as random between one annihilation event and another. Experimental 
observables should then be the average over all the different directions in isospin space. 

To proceed with random grooming, we write the amplitudes h(k) in spherical components 


h aii = J^-Y^(k)H a (k) 





Y lfi {k)H b {k) 


(26) 


with 







(27) 


H a (k) 

H b {k ) 


= ie ik ikG'(k 2 )( 1 + - 
= —ie~ lk ^kG\k 2 )( 1 


v ■ 

v ■ k 

Uk 


Independent grooming of the Skyrmion and anti-Skyrmion means identifying the total am¬ 
plitude with the /i-th component of isospin as 





Vl„{n)Y lv (k)H a (k) + V^(S )Y,„(k)H t (k) 


(28) 


with Q and 5 being two sets of independent Euler angles. Observables calculated from the 
coherent state (|25j) with the amplitude (|28|) are to be averaged over the range of O and 5. 

First we examine the one-body rate, that is the probability for finding a pion of momen¬ 
tum k and isospin type // in the coherent state. We have 


ir ( 1 V,fc) = j dndzh„(k,n,E) 

= —k / dSlrE 
48?r 3 J 

vli(n)vlMH* a (k)H a (k) 

+vi:(z)vi u ,(m:(k)H b (k) 


+KWKAm:(k)H b (k) 


V 


1* / 
[US V 




(29) 


After the angular average, we obtain 


W (1) (ju, k) 


-[H* a (k)H a (k) + H* b (k)H b (k)] 


-k 2 {G\k 2 )f 


, v ■ k 9 

(i + —) 2 + (i 

Uk 


v ■ k 

k ’k 


2 


(30) 


which is independent of /i and the impact parameter b, as we expect of a one-body observable 
averaged over independent groomings. 

Next we look at the 2-body rate, that is the probability of finding a pion of momentum 
k and charge type /i and a pion of momentum q and charge type p in the state \h). We have 
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W (2 \/i,k;p,q) 


j dVdZ /i*(fe,fi,S) hpfan, Z) h* p {q,n,Z) h p {q,n,Z) 

[ dUdZ 
36vr 2 7 

H-^(n)2^^(s)x»^(n)2^^(n)^(fc) J ff 6 (fc) J H-*( flr ) J H a (qr) 

H-^(n)2?^(s)x»^(n)2^^(s)^(fc) J H- 6 (fc)jy ( :( fl r) J H- 5 (qr) 

+Di;(fi)I ^(Z)Vl:(Z)Vl a ,(Q)H:(k)H b (k)H;(q)HM) 
H-^(n)2?^(s)x»^(s)x»^(H)^(fc) J ff 5 (fe)^( fl r)jy 6 (< l ) 

+(a <_> 6, n <- S)] (31) 


Terms with odd numbers of 11 or S average to zero. Applying the following formulae for 
Wigner functions 


87 r 2 

£>!*£>L = — 5 il a 


pis pa 


pp'-'isa 5 


(32) 


D 1 D 1 = (-l)' 1- " / D 1 * D 1 

ill/ pa \ J / — a—u p 


per 


87T 2 


(-Ir^^-A,-*, 


(33) 


fKKK'K' = 8l2 £ 5*1-7(11;PP|AM> 2 (ll; W |AJl/')(H;r'V|AM'>, (34) 

J A MM' ZiV + 1 

the two-body rate is then 

Wr (2) (p, *;*>,«) = i[ A^iH'MH^H’MWq) + (o« 6)) 
+B„(z)(H;(k)H„(k)H;(q)H b (q) + (a « 6)) 
+C w (z)(H;Xk)H b (k)H:(q)H b (q) + (a « 6)) 

+AA) (HHk)H b (k)H- b (q)H b (q) + (o « 6))] (35) 
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where 


A^Z) = 9 £(-l) L (ll; PP\LM) 2 (11; 00|AO ) 2 W(llll; LA)P a (z) 

LA 

= 1 

G'pp(z') p Z 

= ^p,p % 

with z — k ■ q. 

Putting the explicit forms of H a and Hb (P7|) into hh^ 2 \ we have 


( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 


W (2) = i k 2 {G\k 2 )) 2 q 2 (G\q 2 )) 2 [ 

i . .. n ■ k 


, , V ■ k r. 

((1 + - ' 2 

Wk 


^k 

) 2 ( i- 

v ■ k 

Wk 


■) 2 (1 + 

v ■ q 

) 2 

+ (i- 

UJ q 


v q^ 2 

+ (1 

+ 

v ' q ) 2 , 

U q 


UJ q 

-) 2 )d - 

-r 

q 

) 2 )2 COS 


v ' k Y{\ - _: q f) 


C0 n 


^_fc) 2 ) 

Wk 


UJfi 


+ S^pZ 2 (1 - (—-) 2 )(1 - (—-) 2 ) 2 cos((fc — q) ■ b) 

UJk UJ q 

The quantity A pp {z) ([16]) is easily evaluated to be 

~w + if H 7 ^ 0 and p 7 ^ 0 


H f 3 - 2 


— j-z 2 if only fi or p is zero 


+ | z 2 if n = 0 and p = 0 


5 1 5 


We should also observe that 

Ap P (z) = App(z), ^ App(z) = 9 . 


pp 


(40) 


(41) 


(42) 


There are two kinds of pion correlations studied in the literature. The first is the standard 
correlation function defined by 


n , \ w 2 (t*,p;p,q) , 

c 2 (m,p; p,q) = ———- -1 


W 1 ( f i,p)W 1 (p,q 

and the second form, most used in extracting pion correlations from annihilation data, is 


(43) 


11 


















i? 2 (/i,p; v, g/(p,p; o,q) 


W 2 (p,p; g) 
^2(p,p; o-, <?)' 


(44) 


This form is suitable for studying correlations from many random sources, where the two- 
body rates for different type of particles (like 7 r + 7 r - ) approaches the product of the corre¬ 
sponding one-body rates. This form is commonly used for ratios of correlations among like 
to unlike pions (eg n = u). 

Annihilation occurs over a range of the impact parameter b. Substantial annihilation 
is expected for b < lfm while it should be minimal for b much large than lfm. To further 
examine the correlations, we average over the impact parameter b with the probability e~ Xb . 
Here A is a range parameter of the order of one inverse fermi, specifying the significant 
region of collision of Skyrmion and anti-Skyrmion. I 11 the two-body rate m only the last 
two terms depend on b. The integral needed for the average is 



(45) 


where Q = k ± q, Q± is the magnitude of the component of Q in the b plane, and (j) is the 
angle between this component and b. This integral gives 




m 


V. LORENTZ TRANSFORMED HEDGEHOG 

When the velocity of the moving hedgehog is not small compared to the speed of light, 
the Lorentz transformed field configuration should be used. In the overall CM frame, the 


moving hedgehog (in the direction of positive x-axis) (j) a in ([[]) can be written as 
fa = (7(1 -vt),y+ 6/2, z)j 2 g((x - vtf + (y + b/2) 2 + z 2 ) 


(47) 


considering that the pion field transforms as a pseudo-scalar. Here 7 = l/^Jl — v 2 /c 2 . 
Similarly the Lorentz transformed anti-Skyrmion is 
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4>b = -( 7(2 + vt), y - V 2 . z)s(i 2 (x + vt) 2 + (y- V 2 ) 2 + 7 2 ) 


(48) 


Compare (ED, © with the Galilean boosted held configurations [1] and their Fourier 
transforms [TO], we define the following “capitalized” quantities 


V = jv 

R= (7 x,y,z) 

K = (k x /l,k y ,k z ). (49) 


The calculation of the fourier amplitudes for the Lorentz case is essentially the same as in the 
Galilean case, if we replace the lower case vectors in the Galilean case by these “capitalized” 

vectors. We should also note that d 3 r = -d 3 R. 

7 

The Fourier amplitudes for the Skyrmion and antiSkyrmion are then 


-h a (K,u; k ) =i— e^G'(iC 2 )(l + —) 
7 7 u k 

-h a (K,u k ) = -t—e ik -iG\K 2 )(l - — 

7 7 


(50) 

(51) 


Independent random grooming on the (anti)hedgehog is simply a rotation on the vector 
K for the Lorentz case. Consequently, the one-body and two body rates, ([k]) and (|35|), are 
replaced by 
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(52) 


and 


Hf > (/., fc; v, q) = (n,K;v,Q) (53) 

7 

Here L stands for Lorentz and G Galilean. The cosine of the relative angle in (|35p should 
be replaced by K ■ Q. In the plane perpendicular to the velocity, the two cases give identical 
pion correlations. 
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VI. RESULT AND DISCUSSION 


Pion correlations arise from the incoherence in emission sources, commonly from thermal 
or experimental averages. In the present case, however, the Skyrme dynamics (in particular, 
the coupling between spatial and isospin degrees of freedom) also leads to interesting cor¬ 
relations. Although related to the regular type of correlations the present study points 
to correlations not only among pions of like charges but among unlike charges as well. Our 
present consideration does not include full dynamical evolution nor final state interactions 
among the pions. However, we hope important aspects of the Skyrme physics in NN an¬ 
nihilation are elucidated here and the qualitative features will remain in a more thorough 
study. 

As an example of what our picture gives for correlations, we study the R 2 of (H)- This 
is the correlation function most commonly studied in annihilation. What is usually studied 
is not the full R 2 as a function of both momentum variables, but rather R 2 as a function of 
pion pair relative momentum summed over all total momenta. That is we define Q = p — q 
and P = p + q. Then consider R> for fixed 0 summed over all P as well as being averaged 
over the random groomings and impact parameter. Note that after those averages, R 2 is a 
function of the magnitude of Q only. We study R 2 using the Lorentz boosted form for two 
values of the incident energy, or velocity, v = 0.1c and v = 0.9c. In the sum over P for fixed 
Q we take into account the constraints of energy and momentum conservation. In order to 
do this calculation we need to choose a form for G[k 2 ) defined in (15). We take the simple 
form, G{k 2 ) = 1 /{k 2 + p 2 ), and take p = A = lfm -1 , where A is the range in the impact 
parameter average. We study R 2 as a function of Q for various charge ratios. These are 
shown in Figure 2. For all cases R 2 tends to one for large Q, reflecting the absence of charge 
correlations for large relative momentum. The range over which R 2 differs significantly from 
one is lfm -1 and is set by our choice of scale. In all cases the like charge pion to unlike 
charge correlations resemble closely the correlations seen in the data and usually attributed 
to Bose-Einstein effects. Note that we see them here although we have made no statistical 
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or thermal assumptions. As we have discussed before ||, we have averaged over P, impact 
parameter and grooming. 


Similar correlations (especially among unlike charged pions) has been discussed [8j for 
heavy ion collisions as a consequence of Skyrme dynamics, ffowever, it may be easier to 
examine these correlations experimentally in NN annihilation since the dynamics is much 
simpler. In further extending our study of the manifestation of large QCD in annihilation, 
we need to study quantum corrections, time dependent rotational Skyrmions and their col¬ 
lisions, interactions among the emitted pion waves and the construction of a coherent state 
representing a half integer spin soliton. These considerations may hold the clue to two-pion 
charge-spin asymmetry [|6|] and other interesting experimental findings |ll|. These studies 
will extend the domain of the classical/quantum approach to annihilation physics. 


This work is partially supported by the United States National Science Foundation. 
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FIG. 2. The two pion correlation, R -2 as defined in (|44|), as a function of pion relative momen¬ 
tum Q after grooming and averaging over total momentum of the two pions and the direction of 
antiproton velocity (as discussed in the main text). We show here for different charge combi¬ 
nations. Solid curves are for v = 0.1 c and dashed for v = 0.9c. 
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